WavePacket: A Matlab package for numerical quantum dynamics. III. Quantum-classical simulations and surface hopping trajectories

WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I (Schmidt and Lorenz, Comput. Phys. Commun. 2017, 213, 223] and Part II (Schmidt and Hartmann, Comput. Phys. Commun. 2018, 228, 229] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagation techniques to WavePacket. There classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces, trajectories may switch between them. To model these transitions, two classes of stochastic algorithms have been implemented: (1) Tully's fewest switches surface hopping and (2) Landau–Zener-based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring nonadiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.1.0, which is essentially an object-oriented rewrite of previous versions, allowing to perform fully classical, quantum-classical and quantum-mechanical simulations on an equal footing, that is, for the same physical system described by the same WavePacket input. The software package is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics are available. © 2019 Wiley Periodicals, Inc.     

Measuring fluorescence to track a quantum emitter's state: a theory review

We review the continuous monitoring of a qubit through its spontaneous emission, at an introductory level. Contemporary experiments have been able to collect the fluorescence of an artificial atom in a cavity and transmission line, and then make measurements of that emission to obtain diffusive quantum trajectories in the qubit's state. We give a straightforward theoretical overview of such scenarios, using a framework based on Kraus operators derived from a Bayesian update concept; we apply this flexible framework across common types of measurements including photodetection, homodyne, and heterodyne monitoring and illustrate its equivalence to the stochastic master equation formalism throughout. Special emphasis is given to homodyne (phase-sensitive) monitoring of fluorescence. The examples we develop are used to illustrate basic methods in quantum trajectories, but also to introduce some more advanced topics of contemporary interest, including the arrow of time in quantum measurement, and trajectories following optimal measurement records derived from a variational principle. The derivations we perform lead directly from the development of a simple model to an understanding of recent experimental results.     

沙丘背风侧不同铁路结构形式对风沙环境的适应性分析

当铁路穿越大风沙漠地区时,风沙灾害会对铁路工程及其正常运营产生严重威胁,而设计一种合理的铁路结构形式能够减小风沙沉积对铁路工程的危害. 在本文中,以敦煌至格尔木铁路沙山沟段落为研究对象,采用多相流的方法对越过沙丘的风沙运动过程进行数值模拟,分别讨论了风沙运动对位于沙丘背风坡的铁路路基工程和桥梁工程的影响. 主要的模拟结果显示:路基工程明显降低了风速并且将沙丘后的回流区分成了两部分,而桥梁工程的导流效应则压缩了沙丘背风坡的回流区;轨道间的道碴增大了铁路表面的粗糙度,在轨道间有少量沙粒沉积,而路基工程两侧则有大量积沙;铁路表面的积沙量与摩阻风速呈现出非线性关系,随着摩阻风速的增大,路基工程沙粒沉积的增加速度大于风蚀能力的增加速度,而桥梁工程则正好相反. 在防止风沙危害铁路方面,设置桥梁工程明显优于路基工程. 本研究为风沙运动对铁路工程的影响提供了理论支持,也为今后的铁路工程设计提供了新的思路与研究工具.      

垂直旋转流场中单颗粒运动状态判别及受力分析

为研究单颗粒在旋转流场中的运动状态及受力情况,以毫米级球形颗粒为例,利用旋转流场颗粒运动装置,通过使用摄像机记录颗粒在流场中的运动轨迹以获取其运动参数,分析了不同转速和颗粒直径条件下颗粒的运动轨迹,拟合得到了颗粒运动状态判别公式以及颗粒运动轨迹公式,分析了颗粒在旋转流场中的受力情况。结果表明,颗粒在旋转流场平衡状态下运动状态主要分为两类,一类是未离开壁面保持静止,另一类是离开壁面保持稳定周向运动;颗粒进行周向运动的轨迹为椭圆形,并且圆心随着转速的增大靠近旋转中心,而随着粒径的增大靠近壁面;颗粒在旋转流场的运动过程中主要受到离心力和旋转科式力作用。     

Experimental study on age and gender differences in microscopic movement characteristics of students

Campus security has aroused many concerns from the whole society. Stampede is one of the most frequent and influential accidents in campus. Studies on pedestrian dynamics especially focusing on students are essential for campus security, which are helpful to improve facility design and emergency evacuation strategy. In this paper, primary and middle school students were recruited to participate in the single-file experiments. The microscopic movement characteristics,including walking speed, headway, gait characteristics(step length, step frequency and swaying amplitude) and their relations were investigated. Age and gender differences in the headway-speed diagram and space requirements were analyzed by statistical tests. The results indicated that the impacts of age and gender were significant. There were three stages for the influence of gender on the headway-speed diagram for both age groups. The impacts on students' space requirements were consistent for different age and gender groups. But the impacts of age and gender on free-flow speed were affected by each other. Due to the connection of walking speed and gait characteristics, the comparisons of gait characteristics between different ages and genders were performed to understand the corresponding differences in speed more deeply. The results showed that differences in step length and swaying amplitude between males and females were significant for both age groups. The effect of gender on step frequency was significant for primary students. But for middle school students,whether gender had significant impact on step frequency was not clear here because of the large P-value. Besides, the influence of age on gait characteristics changed with gender.     

Scaling of particle trajectories on a lattice

The scaling behavior of the closed trajectories of a moving particle generated by randomly placed rotators or mirrors on a square or triangular lattice is studied numerically. On both lattices, for most concentrations of the scatterers the trajectories close exponentially fast. For special critical concentrations infinitely extended trajectories can occur which exhibit a scaling behavior similar to that of the perimeters of percolation clusters.At criticality, in addition to the two critical exponents =15/7 andd _f=7/4 found before, the critical exponent =3/7 appears. This exponent determines structural scaling properties of closed trajectories of finite size when they approach infinity. New scaling behavior was found for the square lattice partially occupied by rotators, indicating a different universality class than that of percolation clusters.Near criticality, in the critical region, two scaling functions were determined numerically:f(x), related to the trajectory length (S) distributionn _s, andh(x), related to the trajectory sizeR _s (gyration radius) distribution, respectively. The scaling functionf(x) is in most cases found to be a symmetric double Gaussian with the same characteristic size exponent =0.433/7 as at criticality, leading to a stretched exponential dependence ofn _S onS, n_Sexp(–S ^6/7). However, for the rotator model on the partially occupied square lattice an alternative scaling function is found, leading to a new exponent =1.6±0.3 and a superexponential dependence ofn _S onS.h(x) is essentially a constant, which depends on the type of lattice and the concentration of the scatterers. The appearance of the same exponent =3/7 at and near a critical point is discussed.     

Copenhagen quantum mechanics

In our quantum mechanics courses, measurement is usually taught in passing, as an ad-hoc procedure involving the ugly collapse of the wave function. No wonder we search for more satisfying alternatives to the Copenhagen interpretation. But this overlooks the fact that the approach fits very well with modern measurement theory with its notions of the conditioned state and quantum trajectory. In addition, what we know of as the Copenhagen interpretation is a later 1950s development and some of the earlier pioneers like Bohr did not talk of wave function collapse. In fact, if one takes these earlier ideas and mixes them with later insights of decoherence, a much more satisfying version of Copenhagen quantum mechanics emerges, one for which the collapse of the wave function is seen to be a harmless book keeping device. Along the way, we explain why chaotic systems lead to wave functions that spread out quickly on macroscopic scales implying that Schrödinger cat states are the norm rather than curiosities generated in physicists’ laboratories. We then describe how the conditioned state of a quantum system depends crucially on how the system is monitored illustrating this with the example of a decaying atom monitored with a time of arrival photon detector, leading to Bohr’s quantum jumps. On the other hand, other kinds of detection lead to much smoother behaviour, providing yet another example of complementarity. Finally we explain how classical behaviour emerges, including classical mechanics but also thermodynamics.     

Novel permutation measures for image encryption algorithms

This paper proposes two measures for the evaluation of permutation techniques used in image encryption. First, a general mathematical framework for describing the permutation phase used in image encryption is presented. Using this framework, six different permutation techniques, based on chaotic and non-chaotic generators, are described. The two new measures are, then, introduced to evaluate the effectiveness of permutation techniques. These measures are (1) Percentage of Adjacent Pixels Count (PAPC) and (2) Distance Between Adjacent Pixels (DBAP). The proposed measures are used to evaluate and compare the six permutation techniques in different scenarios. The permutation techniques are applied on several standard images and the resulting scrambled images are analyzed. Moreover, the new measures are used to compare the permutation algorithms on different matrix sizes irrespective of the actual parameters used in each algorithm. The analysis results show that the proposed measures are good indicators of the effectiveness of the permutation technique.     

薄板结构加筋布局优化设计方法研究

提出一种两步优化策略,以加筋板结构的固有频率最大化为目标函数,以结构所受外载荷作用的最大静变形为约束条件,开展薄板结构加筋构件的布局优化设计研究。为了降低加筋布局优化的难度,提高优化设计的效率,将加筋等效为一系列弹性铰（点）支撑,以便快速获得加筋横向移动的灵敏度信息。在基本不改变结构重量的情形下,通过合理布局加筋位置,能显著改善结构的刚度分布,提高结构的整体承载能力。随后,小幅调整加筋的截面尺寸,以满足对结构最大变形的设计要求。最后,用两个算例验证了所提优化方法的可行性和有效性。     

An adaptive method for computing invariant manifolds in non-autonomous, three-dimensional dynamical systems

We present a computational method for determining the geometry of a class of three-dimensional invariant manifolds in non-autonomous (aperiodically time-dependent) dynamical systems. The presented approach can be also applied to analyse the geometry of 3D invariant manifolds in three-dimensional, time-dependent fluid flows. The invariance property of such manifolds requires that, at any fixed time, they are given by surfaces in R³. We focus on a class of manifolds whose instantaneous geometry is given by orientable surfaces embedded in R³. The presented technique can be employed, in particular, to compute codimension one (invariant) stable and unstable manifolds of hyperbolic trajectories in 3D non-autonomous dynamical systems which are crucial in the Lagrangian transport analysis. The same approach can also be used to determine evolution of an orientable ‘material surface’ in a fluid flow. These developments represent the first step towards a non-trivial 3D extension of the so-called lobe dynamics — a geometric, invariant-manifold-based framework which has been very successful in the analysis of Lagrangian transport in unsteady, two-dimensional fluid flows. In the developed algorithm, the instantaneous geometry of an invariant manifold is represented by an adaptively evolving triangular mesh with piecewise C² interpolating functions. The method employs an automatic mesh refinement which is coupled with adaptive vertex redistribution. A variant of the advancing front technique is used for remeshing, whenever necessary. Such an approach allows for computationally efficient determination of highly convoluted, evolving geometry of codimension one invariant manifolds in unsteady three-dimensional flows. We show that the developed method is capable of providing detailed information on the evolving Lagrangian flow structure in three dimensions over long periods of time, which is crucial for a meaningful 3D transport analysis.